function diff = model_shock(theta,mean,param,aux)

    param.lambda = exp(theta);
    param = model(param,aux);
    diff=(mean_fun(aux.n,'x',param)-mean);

end